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ABSTRACT 



Context. Previous SiO maps of the innermost regions of HH212 set strong constraints on the structure and origin of this jet. They rule 
out a fast wide-angle wind, and tentatively favor a magneto-centrifugal disk wind launched out to 0.6 AU. 

Aims. We aim to assess the SiO content at the base of the HH212 jet to set an independent constraint on the location of the jet launch 
zone with respect to the dust sublimation radius. 

Methods. We present the first sub-arcsecond (0744 x 0796) CO map of the HH212 jet base, obtained with the IRAM Plateau de Bure 
Interferometer. Combining this with previous SiO(5^4) data, we infer the CO(2-l) opacity and mass-flux in the high-velocity jet and 
arrive at a much tighter lower limit to the SiO abundance than possible from the (optically thick) SiO emission alone. 
Results. Gas-phase SiO at high velocity contains at least 10% of the elemental silicon if the jet is dusty, and at least 40% if the jet is 
dust-free, if CO and SiO have similar excitation temperatures. Such a high SiO content is challenging for current chemical models of 
both dust-free winds and dusty interstellar shocks. 

Conclusions. Updated chemical models (equatorial dust-free winds, highly magnetized dusty shocks) and observations of higher J 
CO lines are required to elucidate the dust content and launch radius of the HH212 high- velocity jet. 
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The jet origin in accreting young stars and its impact on star and 
planet formation is still a major enigma. Although a magneto- 
hydrodynamical (MHD) process appears to be required, it is still 
debated which fraction of the mass-flux originates from the stel- 
lar surface, from reconnections in the magnetosphere-disk inter- 
action zone, or from the disk surface, with all three regions prob- 
ably contributing to some degree (see eg. Ferreira et al. 120061 for 
a review). Rotation searches in atomic microjets from T Tauri 
stars set an upper limit of =* 0.2-3 AU on the launching radius 
of any atomic steady magneto-centrifugal disk wind (Anderson 
et al. 120031 Ray et al. 120071 and references therein). On the other 
hand, MHD disk winds launched beyond 0.2-1 AU are expected 
to be mostly molecular (Panoglou et al. 120121 1. Therefore the 
molecular component of jets holds an important clue to the total 
radial extent of the launch zone. 

A key species in this respect is SiO, a specific tracer of 
molecular jets from the youngest Class protostars, with min- 
imal contamination by ambient swept-up gas (Guilloteau et al. 
[19921 Guefh et al. [19991 Hiran o et al [20061 [20TD1 >. The HH212 
jet in Orion (Zinnecker et al. 1998) provides particularly in- 
teresting constraints on jet structure. An inner bipolar microjet 
peaking at ± 500-1000 AU from the source was discovered in 
SiO (Codella et al. l2007l hereafter Paper I). The fastest SiO gas 
reaches radial velocities +10 km s from systemic, ie. a flow 
speed V p - 150 km s _1 after correction for inclination (^4° from 
the plane of the sky; Claussen et al. 119981) . arguing that it is trac- 
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ing material closely associated with the primary jet. This view 
is supported by its narrow width < 90 AU (Cabrit et al. 2007 
hereafter Paper II) and by the small blue/red overlap in SiO de- 
spite a quasi edge-on view, which restricts the opening angle of 
the fastest SiO gas to < 4°-6° (Paper I). At low radial velocities 
less than ±6 km s _1 from systemic, SiO traces broader struc- 
tures consistent with expanding jet-driven bowshocks (Lee et al. 
2008). The pointed shape of these bowshocks requires that the 
wind speed drops sharply away from the axis. This appears to be 
inconsistent with the fast wide-angle wind predicted by the X- 
wind model (cf. Cai et al. 120081 ). and requires a slow wide-angle 
wind (if any) more compatible with an extended disk wind. 

An extended MHD disk wind also seems favored by the ten- 
tative rotation signatures reported across the tip of the SiO bow- 
shocks. A magneto-centrifugal launch radius of 0.3-0.05 AU 
was inferred assuming a flow speed V p =100-200 km s~'(Lee 
et al. 2008); but the low mean radial velocity of the rotating gas, 
\V— V S ysl = 1 -5-4 km s , implies a lower deprojected flow speed 
V p ^20-60 km s , hence a larger launch radius of 0.6 AU in 
both lobes (see Equ. 5 in Anderson et al. 2003), if the rotation 
interpretation is correct. Unfortunately, no rotation estimates are 
yet available for the fastest SiO jet material, as it remains unre- 
solved laterally at 073 resolution (Paper I; Lee et al.[2008 ). 

The SiO gas-phase abundance may offer an independent clue 
to the location of the jet origin with respect to the dust sublima- 
tion radius (^ 0.2 AU), since silicon is one of the main grain 
constituents. Our previous line ratio analysis in the HH212 mi- 
crojet shows that SiO is optically thick for J up = 2 to 8. The im- 
plied minimum SiO abundance strongly depends on the poorly- 
known jet density and line-of-sight velocity gradient, and varies 
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between 0.05% and 90% of the elemental silicon (Paper II). In 
this letter, we present new PdBI CO(2-l) observations of the 
HH212 microjet, which allow us to put more severe constrains 
on its density, mass-flux, and SiO abundance. We compare the 
latter with chemical models of dust-free vs. dusty winds, and 
outline the work required to proceed in solving this question. 

2. Observations 

CO(2-l) observations of the central region of the HH212 out- 
flow were obtained in February 2008 and January, March, and 
April 2009 with the IRAM Plateau de Bure Interferometer 
(PdBI) in France. The six-element array was used in its A, B, 
and C configurations (baselines from 48 m up to 760 m) for a 
total time of -23 hours. The CO(2-l) line at 230.5380 GHz was 
observed with a 40 MHz (~ 50 km s~') bandwidth and a chan- 
nel sampling of ~ 0.1 km s , later smoothed to 1 km s to 
increase the signal-to-noise ratio. Two units with a bandwidth of 
320 MHz were used to measure the continuum. The data were re- 
duced using the public GILDASQ software. Amplitude and phase 
were calibrated on observations of 0528+134 and 0605-085. 
The absolute flux density scale was determined on MWC 349 
with an estimated uncertainty of ~ 25%. Line and continuum 
maps were produced using natural weighting and were restored 
with a clean beam of 0796 x 0744 (PA=20°). Since we are inter- 
ested only in the innermost jet knots within ±2" of the central 
source, no correction was applied for primary beam attenuation 
(HPBW si 22"). The flux, position, and size of the continuum 
source all agree with the 1 .4mm data presented in Paper I, ex- 
cept for the tentative source MM2, which is not confirmed. The 
secondary source seen at 850/im by Lee et al. (2008 ) at Aa — 
+ 1", AS = -077 is also unconfirmed (see Fig. [TJ. In the follow- 
ing, we concentrate on our results in CO(2-l). Radial velocities 
are expressed in the V L sr rest frame unless otherwise specified. 

3. Results 

3.1. The inner 5": outflow cavity vs. jet emission 

Figure [1] compares in two different velocity ranges the cleaned 
CO(2-l) map of the inner 5" of the outflow with the SiO(5- 
4) map from Paper I, restored here with the same synthe- 
sized beam as for CO(2-l). In the left-hand panel (Fig. [T^) 
both lines are summed over a medium-velocity interval (here- 
after MV) covering [-9.5km s -1 ,— 1.5km s -1 ] in the blue and 
[+4. 5km s~',+9.5km s -1 ] in the red. It can be seen that the CO 
MV emission peaks along the SiO axis, but also delineates a 
broader component not seen in SiO, with a biconical morphol- 
ogy opening away from the exciting source. This corresponds 
to the base of the CO swept-up cavity mapped on larger scale 
by Lee et al. (2006] 120071 1. Owing to this extended compo- 
nent and the lack of short spacings in our u — v coverage, the 
reconstructed CO brightness in the MV range is very sensi- 
tive to the chosen cleaning parameters. In the right-hand panel 
(Fig. [3J>) both lines are summed over a high-velocity interval 
(hereafter HV) covering [-16.5km s _I , -9.5km s _I ] in the blue 
and [+9. 5km s _I ,+12.5km s -1 ] in the red. In this HV interval, 
CO is seen to trace the same gas as SiO, ie. a narrow jet still 
unresolved in the transverse direction, with no significant con- 
tamination from the broad cavity. We verified that contrary to 
the MV emission, the HV emission is compact enough that its 
flux is robust and independent of cleaning parameters. 



HH212: SiO blue peak 




Fig. 2. Top panel: Line profiles in CO(2-l) (solid) and SiO(5- 
4) (dashed) toward the northern SiO peak. Bottom panel: Ratio 
of CO(2-l) to SiO(5-4). The vertical dot-dashed li ne mar ks the 
ambient LSR velocity (+1.6 km s ; Wiseman et al. 2001 ). 



3.2. CO(2-1) opacity in the HV jet 



As shown in Paper II and confirmed by Lee et al. (2008), SiO(5- 
4) is close to LTE and optically thick in the inner HH212 knots. 
On the other hand, the high ratio CO(3 -2)/CO(2-l) > 1.5-2 
seen at high velocity (Lee et al. 2007) requires that CO(2-l) 
is optically thin in the HV jet. We may then infer its opacity 
r^' from the line temperature ratio R of CO(2-l) to SiO(5^4) 
as R =i t1 1 (T^°/T^°), where r e c x ° and 7^° are the excitation 
temperatures of CO and SiO, respectively. Since CO has a lower 
dipole moment than SiO, it is even more easily at LTE and we 
expect that r£° T^. Furthermore, the SiO knots are sufficiently 
far from the protostar (500 AU) and from any other radiation 
source that mechanical heating (eg. shocks, mixing-layers, am- 
bipolar diffusion) will largely dominate over radiative heating 
(Shang et al. 2002 Panoglou et al. 12013 . Therefore, any strong 
local gradient in 7\ should come with a steep accompanying 
velocity gradient, and we will minimize the difference between 
r£° and T^° caused by such gradients by examining the ratio R 
as a function of velocity. 

In Figure [2] we plot the main beam temperature ratio R as a 
function of velocity toward the (northern) blue SiO inner knot. 
In the HV blue range V < -9.5km s 1 , the CO(2-l)/SiO(5- 
4) ratio is characterized by a low value R 0.2, and we infer 
that t1 1 (T^/T^°) 0.2. In the redshifted HV range, the ra- 
tio R appears similar but CO(2-l) is affected by a strong dip 
around Vlsr- +9 km s _1 , also seen in CO(3-2) (Lee et al. l20071 >. 
This dip is present at all positions and may be caused by an ex- 
tended foreground component fully resolved-out by interferome- 
ters. Given this complicating circumstance, we focus on the blue 
HVjet below. 



3.3. HV jet density and mass-flux 



http://www.iram.fr/IRAMFR/GILDAS 



Since the HV blue material is optically thin (Lee et al. 120071 1. 
its beam-averaged CO column density, /Vco(HVB), can be in- 
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Fig. 1. CO(2-l) emission (blue and red contours) overlaid on top of SiO(5-4) emission from Paper I (color scale with white 
contours). Both maps were cleaned with the same beam of (X'96 x 0744 at PA = 20°(cf. filled ellipse in lower left corner), a (left 
panel) Medium-velocity (MV) emission, integrated over [-9.5, -1.5] km s _1 (blue), and [+3.5,-1-9.5] km s _1 (red). CO contours 
range from 6cr (13 K km s ) to 78<xby 6cr. SiO contours range from 6<x (74 K km s _1 ) to 18<x by 3<x. b (right panel) High-velocity 
(HV) emission integrated over [-16.5,-9.5] km s _1 (blue), and [+9.5,+ 12.5] km s _1 (red). CO contours range from 6<x (1.4 and 
1.1 K km s _1 for blue and red, respectively) to 13cr by l<x. SiO contours range from 6<x (21.5 K km s~') to 21cr by 3<x. In both 
panels, 1.4mm continuum from the central source MM1 is shown in light-blue contours, ranging from lOcr (5 mJy beam -1 ) to 70<x 
by steps of lOcr. 



ferred from the CO(2-l) integrated intensity, without opacity 
correction. At the northern peak, we measure T m \,dV = 8.2 K 
km s _I over the HV range and infer that A^co(HVB) ^ 5 x 
10 15 (r e c x o /50/r) cm -2 for LTE in the range ^ 50 - 500 K de- 
termined in Paper II. 

We may also infer a lower limit to the volume density 
of the HV blue jet at the northern peak. With a jet diameter 
D < 072 = 90AU (cf. Paper II), and a beam size across the 
jet b ± = 0744 = 200 AU, we find n H ~ (b±/D 2 )N C o/X C o > 
10 5 (T™/50K) x (10- 4 /Xco)cnr 3 , where X co is the CO gas- 
phase abundance relative to H nuclei. 

Assuming a steady flow along the jet axis, we may 
also estimate the one-sided jet mass-flux as given by 
(see Lee et al. l2007l l M, = lAm H Vjb ± N C oX co - 

10- 7 (r e c x °/50K) (lO- 4 /X co ) M yr-', where Vj = 150 km s" 1 is 
the relevant deprojected speed for the HV jet (cf. Sect. [1). This 
Mj is ten times lower than the value of Lee et al. (120071 1 for the 
same 7^.°, because we excluded the MV CO emission (domi- 
nated by the slower bowshocks and the large-scale cavity). The 
appropriate choice of Xco depends on the assumed dust content 
in the jet. Chemical models of dense Class jets show that CO 
self-screens efficiently against FUV photons and contains most 
of the available gas-phase carbon not locked in grains (Glassgold 



et al. 119911 Panoglou et al. 120121 1. Correspondingly, we expect 
^co - 3.5 x 10 4 (solar carbon abundance; Holweger 2001) 
when the wind is dust- free, and Xco - 10~ 4 for an interstellar 
dust/gas ratio ^ 1% (Flower & Pineau des Forets 2003). Hence, 
the above values for Mj and tin are decreased by a factor ^3.5 
if the HV jet is grain-free, while they both increase by a factor 
10 for a high r e c x ° = 500 K. 

3.4. SiO abundance in the high-velocity jet 

Using the on line RADEX code (Van der Tak et al. [20071 in 
the large velocity gradient (LVG) approximation, we compared 
the LTE opacities of SiO lines versus that of CO(2-l). For the 
SiO(2-l) line and T k = 50-500 K, we find that 

Xs,o = Af(SiO)/AV ^ r 2 s l /T 2 c l /rg°\ 19 
Xco N(CO)/AV (122 ±10) \Tg> ) ' 

With our result from Sect. EH that T 21 (r e c x ° /T^ ) R = 0.2, 
and the observational constraint Tj i 1 > 1 (cf. Paper II), we in- 
fer that Xsio/Xco > 0.04f (r e s x °/r e c x °) - 9 . Adopting the X co val- 
ues discussed above and a solar elemental abundance for silicon 
[Si/H] Q =s 3.5 x 10~ 5 (Holweger. 120011 1. we obtain that gas-phase 
SiO in the inner blue jet knot of HH 212 represents at least 
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70% of the elemental silicon if the wind is dusty, and at least 
=2 40% of it if the wind is dust-free, with an uncertainty factor 
of (T^f/T^®) 9 - The uncertainty in Xsio from this method is ^ 
factor 2 in planar C-shock models. 

4. Implications for the launch radius of the SiO jet 

An SiO abundance reaching > 10%-40% of elemental silicon is 
challenging for current models of MHD winds — both dust-free 
and dusty. We discuss each of these two cases below. 

The chemistry of dust-free MHD winds was investigated 
in ID by Glassgold et al. (1991). The wind is launched from a 
5000 K protostar of radius R - 10R Q = 0.045 AU, and under- 
goes various degrees of acceleration and expansion. The frac- 
tion of silicon in the form of SiO reaches ^ 50% only when 
the base density exceeds a critical value, corresponding to an 
isotropic mass-flux rate > 3 x 10~ 6 M o /yr in their accelerating 
model. This is well above our estimate of Mj < 3 x 10~ 7 M o yr _1 
for the HV jet in HH212, if dust free with T k < 500 K (see 
Sect. 13.3b . Furthermore, when a fiat far-ultraviolet (FUV) ex- 
cess below 2000A is present, SiO is entirely photodissociated 
throughout the wind (see Fig. 10 of Glassgold et al. 1991). The 
assumed FUV photon flux at the SiO dissociation limit ^8eV 
is similar to that produced by an accretion shock of blackbody 
temperature 10,000 K and total luminosity 3.5L . The Class 
protostar of 14L driving the HH212 outflow (Claussen et al. 
1998 ) should have a FUV flux from accretion at least as high, 
hence negligible SiO in the wind. Increasing the launch radius 
would not solve the problem, since both the base wind density 
and the radiation field would drop as l/R^, keeping the same 
ratio of SiO reformation to photodissociation rates. Therefore, 
spherical dust-free winds do not seem able to reproduce the high 
SiO content in the HH212 jet. 

A possible factor favoring SiO synthesis in dust-free winds 
would be if the jet does not arise from but from a narrow 
equatorial annulus of width AR <sc Ro, eg. as assumed in the 
X-wind model (Cai et al. 120081 ). This would increase the base 
density by a large factor Rq/AR for the same jet mass-flux, thus 
enhancing Si + recombination and SiO formation with respect to 
photodissociation. However, this hollow wind geometry would 
also enable the penetration of destructive FUV photons through 
the emptied polar regions (see Panoglou et al. 2012S), an effect 
not included in the ID calculations of Glassgold et al. (I1991I ). 
Detailed chemical wind models in 2D with a FUV accretion ex- 
cess are thus necessary to check whether an equatorial dust-free 
wind could reach an SiO abundance ^ 40% of elemental Si, as 
estimated in the HV jet of HH 212. 

An alternative scenario is that SiO would form by shock pro- 
cessing of silicate grains in a dusty wind launched from the disk 
(Chandler et al. |20%T1 Panoglou et al. l20"T2i ). Recent models of 
C-type shocks with ion-neutral decoupling suggest that sputter- 
ing of charged grains by drifting neutrals can release up to 5% 
of Si in the gas phase, which then reacts with O2 or OH to 
form SiO (Schilke et al. [19951 Gusdorf et al. 2008a). A possi- 
ble problem with sputtering, raised in Paper II, is that the dy- 
namical timescale for the inner SiO knots in HH212 is only 
25 yr, which may be too short to complete the sputtering pro- 
cess and the conversion of Si into SiO (Gusdorf et al. |2008b ). 
Fortunately, at the high densities > 10 5 cm~ 3 of the HH212 in- 
ner knots (cf. Sect. 13.31 ). grain-grain collisions also become im- 
portant (Caselli et al. 119971 ). They will alleviate this problem 
in two ways: First, by shattering dust grains into smaller frag- 
ments, they increase the coupling between charged and neutral 
fluids and shorten the C-shock timescale. Second, they induce 



grain vaporization, which directly releases SiO molecules into 
the gas-phase, reaching up to 5% of the total silicon at 40 km s 
(see Guillet et al. 1201 II ). The same process may also occur in 
single-fluid (J-type) shocks, with SiO reaching 2% of the total 
silicon at 50 km s _1 (Guillet et al. 2009 ). We note that the above 
SiO yields refer only to interstellar shock models with Alfven 
speeds Va of a few km s _1 , a value ^10 times lower than ex- 
pected at 500 AU in protostellar jets (Hartigan et al. 2007 Garcia 
et al. 120011 ). Extrapolating the trends with Va shown in Fig. 5 of 
Caselli et al. ( fT9"9"71 ) and Fig. 6 of Guillet et al. (120091 suggests 
that in dusty MHD jets, grain-grain collisions might be able to 
release 10% of Si in the form of SiO, compatible with our esti- 
mate in HH212 in the dusty wind case. Dedicated shock models 
with appropriate initial conditions would be necessary to verify 
this extrapolation. 

In conclusion, ^ 10%-40% of elemental silicon in the form 
of SiO is challenging for current chemical models of both dust- 
free winds and dusty shocks, and calls for additional modeling 
with more realistic assumptions to clarify the maximum SiO 
abundance that can be reached in either case. ALMA observa- 
tions in CO(3-2) and CO(6-5) will also be essential for better 
constraining the jet temperature, width, mass-flux, and density, 
as well as the ratio T^°/T™ entering our SiO abundance esti- 
mate, all of which are key parameters for distinguishing among 
current MHD jet models. 
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